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Abstract 

It has very recently been demonstrated by Csaki, Kaloper and Terning (CKT) 
that the faintness of supernovae at high redshift can be accommodated by mixing 
of a hght axion with the photon in the presence of an intergalactic magnetic field, 
as opposed to the usual explanation of an accelerating universe by a dark energy 
component. In this paper we analyze further aspects of the CKT mechanism and 
its generalizations. The CKT mechanism also passes various cosmological con- 
straints from the fluctuations of the CMB and the formation of structure at large 
scales, without requiring an accelerating phase in the expansion of the Universe. 
We investigate the statistical significance of current supernovae data for pinning 
down the different components of the cosmological energy-momentum tensor and 
for probing physics beyond the standard model. 



1 Introduction 



The dynamics of the Universe in the standard cosmological model is completely deter- 
mined by the present values Q of the Hubble parameter, Hq, and the ratios f2oi between 
the energy density of the different "matter" components and the critical density. The 
task assigned to cosmologists by Sandage is the precise determination of these numbers. 
Unfortunately, they are not directly measurable and their determination relies on the 
measurement of quantities indirectly dependent on them. Several bounds on the floi 
follow from the study of the anisotropies of the cosmic microwave background radiation 
(CMB), the formation of large scale structure, the age of the Universe, etc. The Hubble 
diagrams (luminosity distances vs. redshifts) of standard candles have traditionally been 
a useful tool for determining constraints on some of the cosmological parameters. One 
of the most astonishing physical results of the past few years is that not only is 
the Hubble diagram of high redshift supernovae (SNe) compatible with previous inde- 
pendent constraints on the Qoi, but it is also complementary to them and in principle 
allows a complete determination {Hq ~ 66 km.s~^.Mpc~^, i^om ~ -3 and i^oA ~ -7). 
In particular this result reveals the presence of a so-called "dark energy" component 
accounting for the dimming of high- 2; SNe as a result of a current period of acceleration 
in the expansion of the Universe. However, both the existence as well as the size of this 
Dark Energy component jeopardize our current understanding of the fundamental laws 
of particle physics and string theory /quantum gravity 0]. 

Very recently Csaki, Kaloper and Terning (CKT) |^ have proposed an interesting 
alternative explanation of combined Type la SNe data from the Calan-Tololo sur- 
vey, the Supernova Cosmology Project (SCP) and the High-2 Supernova Search 
Team (HZT) without requiring a period of accelerating expansion. By assuming 
the existence of a light axion (m^ ~ 10^^^ eV/c^) and an intergalactic magnetic 
field of amplitude |B| ~ 10"^'^ Tesla, the authors of P found that the magnitude- 
redshift curve corresponding to the best fit flat universe with a cosmological constant |@] 
{flom ~ 0.28, flo\ ~ 0.72) can be closely reproduced with cosmic strings replacing the 
cosmological constant. The axion interacts with photons via a term in the effective La- 
grangian of the form, Cint = jjj^ E ■ B, where is the axion field and Ml is the axion 
coupling mass scale, which CKT fit to ~ 4-10^^ GeV/c^. As a result of this coupling, 
photons and axions will oscillate in a background magnetic field via a Primakoff effect. 
Assuming a randomly oriented magnetic field with a coherence length of the order of 
1 Mpc, almost one third of the highest redshift photons would decay by the time they 
reach us. 

In this letter we study the statistical significance of the CKT results, and discuss 
possible generalizations of the CKT mechanism. In order to understand the model 
dependence of generalized photon loss mechanisms we use the SNe data to identify the 
allowed regions in the parameter space of such models, obtaining fits at least as good as 
the cosmological constant accelerating universe. The anisotropy of the CMB still requires 

'^"'^^A subsript will always denote the value taken today by any quantity. 
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a Dark Energy component to the energy-momentum tensor. The large scale structure 
of the Universe then imposes constraints on the equation of state of the invisible Dark 
Energy. The difference from previous analyses is that the SNe data now allows for a 
significantly larger range of the equation of state parameter wx (in p = wxp) for the 
Dark Energy component. It is interesting that the preferred region for wx remains in the 
range wx < 0, in agreement with requirements of large scale structure. The key point 
is that instead of indicating an accelerating expansion of the Universe, the dimming of 
SNe at high redshift would reveal the presence of a light axion which mixes with photons 
in the intergalactic medium. 

Because this subject unites several branches of theory and experiment which some 
readers may not be familiar with, in Section ^ we briefly review the relevant experiments 
and current bounds on the relevant cosmological and particle physics parameters. In 
Section ^ we review the theoretical basis of the CKT photon loss mechanism and its 
generalizations. This mechanism is intrinsically quantum mechanical, as opposed to 
the essentially classical models of photon absorption by dust 0. We then describe in 
Section ^ the constraints placed on the cosmological parameters by the Hubble diagrams 
of Type la SNe and how they are modified in the presence of a photon loss mechanism. 

In Section ^ we perform a statistical analysis of the data and our fits, and demonstrate 
that the CKT model with a best fit photon decay parameter (with f2om = -3, ^lox = 
.7,wx = —1/3) is as good as the f2om = -28, Qqa = -72 cosmological constant accel- 
erating universe scenario. We present some exclusion regions of the parameter space, 
but there is a large range of allowed parameters for a generalized CKT mechanism to 
accommodate the data. 

We conclude with proposals for future experiments to better pin down the Sandage 
numbers and ferret out the nature of the Dark Energy, and perhaps reveal the existence 
of an axion coupled to the particles of the Standard Model or other new physics beyond 
the standard model. 

2 Experiments and bounds 

In this section we summarize experimental bounds on and estimates of the relevant 
cosmological, astrophysical and particle physics parameters. 

2.1 Cosmological Constraints 

Although the energy density ratios Qot are not directly measurable, they affect the time 
evolution of the scale factor of the Universe which in turn controls at least three types 
of measurable quantities: (i) the growth of density perturbations responsible for the 
formation of structure at large scales; (ii) the age of the Universe; (iii) the propagation 
of light between its emission by some known astrophysical sources and its detection 
today. 
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The sequence of structure formation starting from smaller fluctuations and ending 
with clusters of galaxies is well described by cold dark matter models, and the structure 
observed today could well have evolved from the density perturbations revealed by the 
CMB anisotropy at 2; ~ 1100. The presence of Dark Energy should not inhibit the 
growth of density perturbations, which implies that the Dark Energy should have started 
to dominate only recently. This requires an upper bound on wx = Px/px at least as 
strong as < = [|,|]. 

The fluctuations in the CMB can also be used to predict the distribution of matter, 
baryon and Dark Energy in the Universe. The position of the first peak in the power spec- 
trum as a function of multipole moment depends on the total density (fiom + ^ox) and 
also less sensitively on f2om and fiobaryon- The latest determination by the BOOMERANG 
experiment is about f2o = 1-0 ± 0.1 and f2om ~ -33 ± .03 (and only a small fraction of 
baryons) ||T^. The fiat universe is heavily favored by this observation, as any deviation 
from flatness tends to grow quickly in Big Bang cosmology. This observational result as 
well as a nearly scale-invariant spectrum of primeval density perturbations is in accor- 
dance with the predictions of inflationary models. The estimate of fiom seems also in 
agreement, at least for its order of magnitude, with the measurements of mass-to-light 
ratios in clusters. Even if the CMB cannot directly teach us about the Dark Energy (the 
CMB is an image of the Universe at 2 ~ 1000 where the Dark Energy was completely 
washed out), it still points out a missing component of the energy- momentum tensor. 

Some other cosmological constraints, for example from the age of the Universe or 
gravitational lensing, are less predictive due to high uncertainties and limited data. 

Finally, we can determine some cosmological constraints from the Hubble diagrams 
of some "standard candle" sources: if we know the absolute luminosity of the source, 
by measuring the flux we receive today, we can deduce the distance travelled by the 
light, and this distance depends on the VLqi ratios. At very low redshift {z < .1), this 
dependence is universal while at large redshift [z ~ 5) all the components with < 
have been washed out by the matter energy density. The intermediate interval is thus 
the best suited to bring some information on the Dark Energy. Among the astrophysical 
objects observed in this redshift interval are the Type la SNe. And by chance. Type la 
SNe have been found to be decent standard candles by means of matching luminosity 
vs. time curves to the phenomenological Phillips curve . As we will explain in detail 
later, the information on the cosmological parameters inferred from the Hubble diagrams 
relies on the nature of the propagation of light. Assuming a free propagation of the 
photons in the cosmological background, and assuming a flat universe, the best fit of 
the SNe data gives |0] fiom = 0.28 ± 0.09. This is in good agreement with several other 
determinations, which is perhaps why the relatively poor current statistics in high- 2; SNe 
data is sometimes overlooked. The aim of this work is to study how this conclusion is 
modified when the photons are no longer freely propagating but rather interact with an 
axion in such a way that some photons are lost between their emission by the SNe and 
their detection today. The following section will review the CKT photon loss mechanism 
and we will analyze the SNe data taking into account such a photon loss mechanism in 
Section ^. 
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2.2 Astro-particle Constraints 



The axion is a proposed ultralight pseudoscalar field that would couple weakly with 
ordinary stable matter. The relevant coupling for the present study is that between 
axions and photons. The coupling is through a term in the Lagrangian of the form, 
a D ]y^72 E ■ B, where is the axion and Mx, is the axion coupling mass scale. 
Several searches for axions have turned up negative, placing constraints on the axion 
mass and couplings. In particular, a study of conversion of solar axions to X-rays in 
a strong magnetic field give > 1.7 x 10^ GeV/c^ |13[. Based on mixing between 
axions and photons in an external magnetic field, precisely the effect needed for the 
CKT mechanism, studies of globular-cluster stars provide a rough limit Mi^ > 1.7 x 
10^0 GeV/c2 
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For more complete reviews, see ||T5[. Bounds on the axion mass are 
closely tied to the axion coupling, but the absence of intergalactic line emissions and 
other experiments place a rough bound of < 10~^ eV/c^ ||16|. The CKT axion of 
mass ~ 10^^^ eV/c^ and coupling of Ml ~ 4 ■ 10^^ GeV/c^ is within current bounds. 
Some further studies and bounds can be found in |T^. 

The mixing between the axion and the photon requires the existence of the magnetic 
field. Magnetic fields appear to permeate the observable universe over a large range 
of scales, from relatively strong fields in the interstellar galactic medium, to weaker 
fields in the outer envelopes of X-ray clusters. The intergalactic magnetic field can be 
deduced by observations of radiation from relativistic electrons around X-ray clusters. 
These observations suggest an intergalactic magnetic field of around B ~ 10~^^ Tesla. 
(See [1^,19] for recent summaries.) The magnitude of the intergalactic magnetic field can 
also be estimated by other means, including studies of radio emission of distant quasars 
which, assuming a coherence length for the magnetic field of about 1 Mpc, yields roughly 
the same estimate ||2^. The source of these magnetic fields is not well understood, 
but there have been several suggestions made, from inhomogeneous cosmological lepton 
number |21] to expanding quasar outfiows [18|. Based on these results, CKT assumed 
a bulk magnetic field of 10~^^ Tesla, and assumed a coherence length of 1 Mpc. The 
distribution and variation of the magnetic fields in the intergalactic medium is still not 
well understood. For a recent review of the relevant experiments and the significance of 



magnetic fields in the early universe, and for a more complete list of references, see |22 . 
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3 CKT photon loss mechanism 



Following the pioneering works ||23| of Sikivie, the CKT model considers an axion-photon 



coupling in the inter-galactic magnetic field 

Cint = ^^(pE-B (3.1) 

where the scale characterizes the strength of the axion-photon interactions. In a 
region where the magnetic field is approximately constant, the polarization of the photon 
whose electric field is parallel to the magnetic field mixes with the axion to form an 
oscillating system, much like that of the massive neutrino system. The mass eigenstates 
in the presence of the background magnetic field are not the same as the mass eigenstates 
in vacuum, and thus the photons oscillate into axions and vice et versa. A quantum 
mechanical computation gives the probability that a photon has not yet oscillated into 
an axion over a distance I in a uniform magnetic field as ^ , 

A 2 2 

P,_,(/) = 1 - , /J. ,, sin^ (iV^^^- V^^^)l/{2hc)) , (3.2) 

where is the mass of the axion, e is the energy of the photon, /i = |B|/ {M^ c^) and 

X± = {m^^c^ ± y^m^"c^~+4/i^) /2 are the squared- mass eigenvalues. Note that despite 

the fact that A_ < 0, that mode has positive energy when the energy e is taken into 
account. 

If /i <^ m^c^ then there are two distinct energy regimes: 

• High energy photons (e ^ m'^c'^/n), for which the oscillation is maximal and 
achromatic: 

P^_^{1) = 1 - sin^ (/i l/{2hc)) (3.3) 



• Low energy photons (m^c^ <^ e <^ m'^c^/fj,), for which the oscillation is small and 
energy-dependent : 

P,^,(/) = 1-44^ sin^ {myi/{Ahe)) (3.4) 

In the Universe, the intergalactic magnetic field is not uniform but is assumed by 
CKT to vary in domains of size -Z^dom ~ IMpc, which is much smaller than the distance 
to high redshift SNe (z ~ 1 is equivalent to distances about 3 x 10^ Mpc) . CKT assume 

Throughout the paper, we wiU use electromagnetic units such that both the electric and magnetic 
fields have the dimension of squared energy while the axion scalar field has the dimension of energy. 
This means that the magnetic field in SI units is given by Bsi — \B\ / \/4TTeoWc^ . The speed of light 
will be denoted by c, Newton's constant hy Qn, and the vacuum permittivity by eg- 
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a uniform distribution of randomly oriented magnetic fields of constant magnitude over 
the 1 Mpc domains^. In that case high energy unpolarized photons (as defined above) 
approach an equilibrium population evenly divided between the two photon polarizations 
and the axion at large distances. CKT found that the approach to equilibrium is through 
an exponential damping of the number of high energy photons 0]: 

— ^7(0 ~ 3 ~^ 3 ^ ^ dec ^ (3.5) 

where the decay length is approximatively given by: 

8 c« Ml 

Before studying the implications of this photon loss mechanism on the SNe data, it 
is worth mentioning how to evade any unwanted effects on the CMB anisotropy. While 
the photons of the CMB were emitted at energies of the atomic sizes, when the magnetic 
field appears they have already redshifted down to much lower energies about 10"^ eV. 
In this low energy regime, the oscillation mixing is energy suppressed and requiring 
1 — P^_»^ to be less than the CMB anisotropy (10"^) gives a lower bound on the axion 
mass around > lO"^*^ eV §. 



4 Cosmography: luminosity distance vs. redshift 
with and without photon loss 

Let us now study the consequences of this photon loss mechanism on the measurement 
of luminosity distances. A nice review of the cosmological parameters and their deter- 



mination for non-specialists can be found in |]26| . 

A homogeneous and isotropic space, as our Universe appears to be at large scales, is 
described by the Robertson- Walker metric 

ds' = -c'de + Rla\t) ( -^4-^ + r'de' + sin' OdA (4.7) 

\ 1 — kr^ J 

where a{t) is the dimensionless scale factor (by convention a(today) = 1), Ro char- 
acterizes the absolute size of the Universe, while with an appropriate rescaling of the 
coordinates the curvature parameter k takes values +1,0,-1 for closed, flat and open 
universes. 



'^'^'As mentioned earlier, these magnetic fields are poorly understood and may be distributed much 
more irregularly than is assumed here. We thank Hui Li for a discussion of this point. 

It has been noted that the photon loss probability and the above conclusions may be modified in 
the presence of the intergalactic plasma , depending on the plasma density. It was suggested in [ p5[ 
that a small free electron density can in fact weaken the constraints on the axion parameters. In any 
case, we have not included such modifications in this analysis. 
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For energy-momentum sources made of a superposition of perfect fluids specified by 
their energy densities and their pressures Pi, the Einstein equations reduce to the 
Friedmann equation (a dot stands for a derivative with respect to the time coordinate): 

H=^Y.P^-^. where i7 = - (4.8) 

together with the conservation equation: 

Pi + 'i{pi+p^)H = Q (4.9) 

It is common to introduce the ratio of the energy densities to the critical energy density 
for a flat universe: 

Note that J2i ^i + ^curv. = 1 by (|4.8| ). In terms of the redshift parameter z = 1/a — 1 (by 
convention, since oq = I, zq = 0) and assuming a constant equation of state pi = WiPi 
for each component of the stress-energy tensor, the conservation equation ( [4.9|) can be 
integrated with respect to time, giving a relation between the densities pi, the scale 
factor a and the equation of state Wi. This gives, 

p, = po^{l + zf^'+^^\ (4.11) 
The Friedmann equation ( |4.8|) can then be rewritten |^ 



H\z)=Hl{l + zf[l + Y,^^{i^+zY'''''"'-m and 1- 



kc^ 



mi 



(4.12) 



A crucial test of cosmological models comes from Hubble diagrams for standard 
candles by plotting their luminosity distances versus their redshifts. The luminosity 
distance is defined by considering a source located at the comoving distance r^, which 
has emitted some light at a time ts with an absolute luminosity C (energy per time 
produced by the source, normalized to that of a standard candle). This light is detected 
today by a detector located at r = which measures a flux (energy per time per area) 
and a redshift z^ . Taking into account the energy loss due to the redshift of the light and 
the increase of time intervals between the emission and the reception, the conservation 
of energy ensures that, 

C = ATirlRl{l + Zsfj^. (4.13) 



(5)pj.oni now we will omit the subscript on the energy density ratios £7^ evaluating those quantities 
only at the present time. 
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The luminosity distance between the source and the detector is defined by, 

L = ^Tid\T. (4.14) 

The quantities characterizing the source are not independent but can be related to each 
other through the ratios f2i governing the dynamical evolution of the Universe. This is 
why the independent measurement of two SNe observables can tell us about Sandage's 
numbers. The trajectory of light is given by the null geodesic equation, ds^ = 0, which 
when integrated becomes, 

dr r'""^^^ cdt cdz , , 

(4.15) 



Using ([4.12|) this integral equation can be solved to relate the comoving distance of 
the source to its redshift]^ 

rAz,) = Sinn | — ^ / = | (4.16) 

\r.hj, (i + z)vi + >:^a((i+z)i+3-._i)y 

This relation allows one to predict the luminosity distance of the source in terms of its 
redshift: 

dL = {I + Zs) rs{z,) Ro. (4.17) 

So far this expression has been derived assuming that once emitted the light propagates 
freely in all spatial directions. The CKT proposal that photons can disappear into 
axions in the intergalactic magnetic field changes the relation ( [4.17|) for dilzg). Imagine 
indeed that over the distance VsRo a fraction (1 — Pj^^{rsRo)) of the number of photons 
has decayed into axions, so that only a fraction P^^y{rsRo) of the energy emitted by 
the source will be detected and the conservation equation ( |4.13| ) has to be modified 
accordingly: 

C = Anrl Rl (1 + z,f J^/P^^^^rM (4.18) 
Then the expression for the luminosity distance as a function of the redshift becomes: 

(1 + Zs) rsizs) Ro 
dL = ^ ^2 (4-19) 

P-y—fji^^s^^s) Ro) 

where rs{zs) is still given by ( |4.16| ). Usually, the data are not presented in terms of 
diizg) but rather in terms of the magnitude defined as 

m = mo + 5 \og^o{HodL/c) (4.20) 



(^■'The definition of the function 'sinn' depends on the topology of the Universe: sinnx = sinh x, x, sin x 
for an open, flat, closed Universe, respectively. See for instance p%. 
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where mo is the magnitude associated to the absolute luminosity of the source (after 
corrections deduced from the fitting to the Phillips curve, mo is the common quantity 
to all SNe and it is fit to the data). 

We would like to comment briefly on the effect of photon loss on the measured flux. 
The luminosities of distant SNe fluctuate too much in general for such SNe to make good 
standard candles. However, it has been noticed that the absolute luminosities of Type 
la SNe can be determined by fitting the time evolution of the measured SNe luminosity 
to the phenomenological Phillips curve [11]. Roughly speaking, the broader the light 
curves of a Type la SN, the brighter it is intrinsically. This relation is specified phe- 
nomenologically by a stretching factor relating the width of the light curve to its intrinsic 
luminosity. Given the absolute luminosity, together with the measured luminosity, one 
determines a normalized flux T . The effect of photon loss is to reduce the measured 
flux, but it is important that the Phillips curve not change as a result. The first point 
is that the Phillips curve measures brightness on a logarithmic scale as a fuction of the 
age of the SN. Hence, if an energy-independent fraction of photons is lost by a partic- 
ular source, the light curve will simply be shifted vertically. This is the effect we have 
been discussing. The only reason for concern is then whether or not the broadness of 
the observed light curve should change as a result of photon loss, as that is how the 
absolute luminosity of a Type la SN is determined. The luminosity signal is generally 
high enough above background that we believe this will be a negligeable effect. A more 
rigorous analysis should be done, however. 

The CKT model described above provides a physical example of a photon loss mech- 
anism in which the P^^^ probability is directly related to particle physics parameters 
that can therefore be constrained by the SNe data. 

What the SNe data indicate (Figure p is that SNe with z ^ 1 are dimmer than 
can be explained by a universe with only ordinary matter in it {wm = 0), assuming no 
loss of photons. The universe undergoes a period of acceleration if the Dark Energy 
has an equation of state wx < —1/3 (the condition for having an accelerating period 
of acceleration today is wx < —1/(3 — SQm) and if —1/(3 — 3^2^) < wx < —1/3 the 
universe will accelerate in the future). Suppose we choose to assume that there is no 
period of acceleration. Then we can study the general requirements of a photon loss 
mechanism to explain the SNe data. We think it is worth studying such generalizations 
because the authors of make many assumptions regarding observational astrophysics 
which are not well understood. For example, as mentioned earlier, the structure of the 
intergalactic magnetic field may be significantly different than the CKT assumption of 
1 Mpc domains of constant magnitude. Furthermore, there may be other acceptable 
photon loss mechanisms, either via new particle physics or via uniform absorption by an 
unknown intergalactic medium. 

For 2; < 1 we require that the photon loss be so as to force a positive slope for the 
magnitude-redshift curve. For example, in the absence of a photon loss mechanism, the 
dotted curve in Figure ^ must be forced upward. The SNe magnitude is proportional to 
the logarithm of its luminosity, so by ( |4.18|) this magnitude is shifted by the addition of a 
term proportional to log(P-y_^^). If we assume that for low z, P^_»^ ~ l — 6P^^^/{z) where 
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5P^^^{z) ^ 1, then we require that 5P^^-y{z) have a hnear term which can control the 
slope of the magnitude-redshift curve. Otherwise a fit will likely become more difficult. 
The other requirement is that if we want to match onto the single Type la observation 



at 2; ~ 1.7 p8[, then the magnitude-redshift curve must be allowed to slope downward 



by then. The easiest way to accomplish this is to give log(P^^^(z)) a small negative 
slope and positive second derivative near 2; ~ 1. These requirements are satisfied by the 
CKT mechanism. We have also checked that in a variety of phenomenological photon 
loss models (without theoretical basis) for which these two requirements are satisfied, 
it is rather easy to fit the SNe data at least as well as the accelerating universe by a 
cosmological constant. 



5 Statistical fits 

In this section we describe the statistics behind our fits for the parameters in the 
generalized CKT model. We have included a total of 18 SNe from the Calan-Tololo 
Supernova Survey 0] and 36 from the Supernova Cosmology Project [0. Following P 
we removed six outlying data points of their 42 total, corresponding to their Fit C. 



Furthermore, in the fits presented here we have tentatively dropped the z ~ 1.7 point ||28 
because of its large error bars compared to the remaining data. We have checked that 
the results are nearly unchanged when that point is included. 

In our fits we have assumed a two component universe with a fractional density fi^ 
of ordinary matter and VLx of Dark Energy with equation of state px = wxPx- We 
assume the CKT model for photon loss parametrized by the photon decay length Ldec- 
Our theoretical models thus involve five parameters: f2m, f^x, wx, mo and Ldec- So 
far all the statistical analyses of the SNe data have arbitrarly reduced this parameter 
space by choosing L^ec = 00 a priori. As we will see now, the introduction of this extra 
parameter to fit the data changes qualitatively and quantitatively the confidence levels 
in the determination of the cosmological parameters. 

Because the absolute magnitude, mg, is related to the calibration of the measure- 
ments, we have chosen to systematically fit this variable to the SNe data (both low 
and high redshift) with a uniform prior. Hence the number of parameters in the fits 
presented here is always one greater than the effective theoretical parameter space. This 
procedure is slightly different than that which is sometimes performed, in which the 
additional annoyance parameters are integrated over. Because the distribution in mo is 
quite highly peaked, the two procedures are nearly equivalent. There will in addition 
be a small correction to these fits from the error in the stretching factor of the Phillips 
curve (we used the effective magnitudes as presented by the SCP collaboration, which 
already involved fitting the stretching factor). In the various fits presented in this letter, 
some further constraints have been placed on the remaining parameters of the model in 
order to make use of previous experimental results and to limit the dimensionality of 
our confidence plots. 

The likelihood of a model to explain experimental data is measured by a together 
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with the number of data points and model parameters, which are then translated into a 
confidence level about the best fit: 



= - ^TW^^T - <)■ (5-21) 



The correlated error matrix aij can be found at [http : //www-supernova . Ibl . gov| for 
the SCP data . 

The data are displayed together in Figure |l], together with magnitude-redshift curves 
of various models including the CKT model. The data are displayed in terms of residual 
magnitudes with respect to the SCP fiat accelerated universe. The CKT model differs 
only slightly from the SCP model for small z. Note that for 1 < 2; < 5 the CKT 
magnitude dips below the SCP accelerated universe. This fact may allow the models 
to be distinguished if better high- 2; data is obtained: a precision of .05 in magnitude 
at z ~ 1 would allow this distinction between the two scenarios to be made. The 
deviation aX z > 100 should not be taken seriously, as the universal magnetic field was 
likely insignificant before large scale structure formation (but see for example ||2^ for an 
alternative point of view). 

Figure ^ demonstrates the maximum likelihood regions for a fiat universe with the 
CMB preferred value VLm = 0.3. In this first statistical analysis, the effective theoretical 
parameter space is two dimensional only, being spanned by the photon decay length, 
Ldec, and the equation of state of the Dark Energy component, wx- Neglecting a photon- 
axion coupling would have restricted ourselves to the line L^cc = 00 m. the plot, with 
the conclusion that at 99% confidence level the Universe is accelerating. Clearly a non- 
vanishing photon-axion coupling opens up a valley that extends beyond wx < —1/3 and 
thus allows the SNe data to accomodate a Universe which does not accelerate, either 
today or in the future. Also note that even with the enlarged parameter space due to 
the photon loss mechanism, negative values for wx are strongly favored, which is an 
independent result consistent with the constraints coming from the CMB anisotropies 
and structure formation. 

It is worth studying what can be learned from current SNe data while relaxing some 
constraints coming from the other cosmological experiments. Of course we expect a 
larger acceptable region of parameter space with the additional photon loss parameter, 
but in particular it is interesting to determine the most likely region in the enlarged 
parameter space. First, let us vary the matter density ratio, fim, keeping a fiat Universe 
{Vtx = 1 — ^m)- The parameter space is now three dimensional. To arrive at the 2D 
confidence levels plotted in Figure 0, for each point {wx, ^^m), we have minimized with 
respect to Ldec (a minimization over the range of values for Ldcc allowed by the different 
constraints on the photon-axion coupling and the intergalactic magnetic field). The two 
figures illustrate the bias introduced in the statistical analysis when removing the decay 
length as a free parameter. Notice again that while the usual interpretation of the SNe 
data (in the absence of photon loss) prefers a current period of acceleration at the 90% 
level, the plot on the right shows that the CKT mechanism can easily accommodate a 
cosmological evolution without acceleration. 



11 



In the last fit reproduced in Figure ^, we have assumed that the Dark Energy is 
formed of a vacuum energy {wx = ~1) and we have examined the bias in the determi- 
nation of the remaining cosmological parameters. Again, for each point flm) in the 
plot on the right of Figure ^, is first minimized with respect to Ldec- Of course, since 
the dimming of SNe is now a combined effect of acceleration and the CKT photon loss 
mechanism, the confidence level regions are enlarged compared to the usual SNe data 
analysis. 



6 Conclusions 

We have presented an analysis of the SNe data taking into account a photon loss mech- 
anism as recently proposed by CKT. This mechanism offers a new parameter, namely 
the photon decay length, to accomodate the data. The main conclusion of our analysis, 
as already announced by CKT, is that the dimming of SNe at large redshift no longer 
requires an accelerating expansion of the Universe, and the photon-axion coupling opens 
up the window for the equation of state of the Dark Energy towards less negative values. 
We argue that the accelerating universe and the CKT models may be distinguishable if 
better high-z data becomes available. In this way future experiments like SNAP may 
reveal new physics beyond the standard models of both particle physics and cosmology. 
In the new allowed region of the parameter space, the SNe data could not only confirm 
the existence of a Dark Energy component but may also reveal the existence of a light 
pseudoscalar that has so far evaded any direct particle physics detection. 
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Figure 1: Hubble diagram for the Type la SNe The residual magnitudes from 

the best-fit flat cosmology with a cosmological constant and no photon loss (the "SCP 
accelerated universe") are plotted. The horizontal dashed line is the theoretical curve 
for the best-fit fiat accelerated cosmology with no axion-photon mixing (fim, ^x-.'Wx) = 
(.28, .72, —1); the solid line is the theoretical curve for the fiat cosmology with a choice of 
photon-axion mixing parameter (fi^, w^x, -^dcc-f^o/c) = (.3, .7, —1/3, 1/3); the lower 
dashed line shows the effect of turning off the photon-axion coupling with the otherwise 
identical choice of parameters. Finally, the upper dot-dot-dashed line is the remaining 
photon intensity probability on a distance rs{z)Ro, with the probability scale given to 
the right of the plot. The lower left plot emphasizes the very high z behavior; and the 
lower right plot enlarges the region near the farthest observed Type la SN. 
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Figure 2: Prom black to white: 68%, 90%, 95% and 99% confidence regions in the 
(wx-L^ec) plane for a flat Universe formed of matter with = .3 and an additional 
Dark Energy component with an equation of state wx- For a given amplitude of the 
intergalactic magnetic field and a given length of coherence for it, the decay length 
is proportional to the square of the axion coupling scale (the values reported on 
the right correspond to B — 10"^^ Tesla and Ldom = 1 Mpc). The dashed and dot- 
dashed lines separate regions of present acceleration (3), present deceleration but future 
acceleration (2) and eternal deceleration (1). 



16 



.2 .4 .6 .8 1 



.2 .4 .6 .8 1 



Figure 3: Best-fit, from black to white: 68%, 90%, 95% and 99% confidence regions in 
the {flm^wx) plane for a flat Universe. The plot on the left is obtained in the standard 
model without mixing between photons and axions; the plot on the right shows the 
modifications due to the photon loss mechanism. The dashed and dot-dashed lines 
separate regions of present acceleration (3), present deceleration but future acceleration 
(2) and eternal deceleration (1). 
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Figure 4: Best-fit, from black to white: 68%, 90%, 95% and 99% confidence regions in the 
(Qm^Q^) plane for a Universe formed of matter and a additional cosmological constant 
The comparision of the two plots shows the influence of the photon loss mechanism (on 
the left the decay length of the photon is set to infinity while on the right it is fit to its 
best value at each point). The dot-dashed line locates fiat universes. 
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